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Abstract 

With a sample of 88.8 x 10^ BB pairs produced at PEP-II in e^e^ annihilation through the 
T{4S) resonance and recorded with the BABAR detector, we search for the B+ — > K*^p^ 
decay mode. A signal is observed for the first time with a significance of more than 5o". 
We measure a preliminary branching fraction of B{B^ K*^p^) = [17.0 ± 2.9(stat) ± 
2.0(syst)l?;9(non-resonant)]xl0-^ The ' 'non-resonant" error corresponds to the uncertainty 
from non- resonant backgrounds not modeled in the fit. The measurement of the longitudinal- 
polarized component to this vector- vector penguin decay is of special interest. We measure 
/i = 0.79±0.08(stat) ±0.04(syst) ±0.02(non-resonant). We measure the direct-CP-violating 
charge asymmetry in this mode to be Acp = —0.14 ± 0.17(stat) ± 0.04(syst). 
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1 INTRODUCTION 



The purpose of this analysis is the simultaneous measurement of the branching fraction, the 
longitudinal-polarization component, and the direct-CP-violating charge asymmetry for the 
B vector- vector decay mode — ^ K*^ p^, which has previously not been seen experimen- 
tally. The measurement of the polarization of this penguin decay is of special interest. There 
is no tree contribution to the decay amplitude and only a very small annihilation process. 
The polarization of the K*^p^ decay mode can be compared to that already mea- 

sured for the similar penguin decay B K*(j). The direct-CP-violating charge asymmetry 
is defined as: 

_ N{B~ ^ if*V) - N{B+ K*\+) 

•^CP = =;^n ■ 

N{B- K*^p-) + N{B+ K p+) 

In this penguin decay, a value significantly different from zero could be a hint for new physics. 



1.1 Physics motivation 

In the Standard Model, the study of the penguin decay B^ K*^p^, combined with 
information from other charmless hadronic B —>■ vector-vector decays into K*p and pp, 
allows us to constrain the angles a and 7 of the unitarity triangle, in a way similar to the 
study of the B —>■ iiii and B —>■ Kn decays pP . The pp modes are used to constrain the angle 
a. The methods that constrain the angle 7 are based on isospin symmetry relationships 
relating the amplitudes of the different K*p decay modes, as well as on the SU{3) flavour 
symmetry relating the K*p and pp modes. They also rely on the relationship a = vr — /? — 7 
and the measured value of the angle f3. A constraint on 7 can already be obtained using only 
the charged-S decays B+ K*°p+, B+ K*+p°, and B+ p+p° 0. The two modes 
B~^ K*^p^ and B~^ P^P^ have been studied previously |Sl E] . A stronger constraint on 
7 can be obtained by including the B^ K*~^p~ mode, a first study of which is also being 
presented at this conference and is based on a similar analysis technique 

Measurements of the rates, polarizations, and direct-CP-violating asymmetries of these 
vector-vector decays permit testing theoretical predictions from the naive factorization jH] 
and QCD factorization [71 |Hj based models. 

In charmless decays of B mesons into two light vector mesons, both longitudinal and trans- 
verse polarization states are possible, but a large longitudinal-polarization fraction, of order 

1 — 4 X ^ 0.9 for K*p, is expected from theory for both tree and penguin decays. Existing 

nig 

measurements of tree dominated vector- vector charmless modes p^p~ and p'^p^ show that 
their longitudinal component is indeed dominant [3 |3J IH]. However, experimental results 
for the polarization in the pure gluonic penguin B^ (j)K*^ and B~^ (pK*^ charmless 
processes indicate a transverse polarization fraction of about 0.5 fOlEllSl- This difference 
can be accounted for in the Standard Model by increasing the non-factorizable contribution 
of annihilation diagrams, through the tuning of some poorly known non-perturbative QCD 
parameters jHl El • The uncertainty on these parameters makes the accurate a priori 
prediction of the polarization difficult for a given decay mode. However, using SU(3) flavour 
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symmetry arguments, the same polarization is expected for the K*^p~^ decay and for 

B (f)K* decays. 

Thus, the comparison of the polarizations obtained in the two penguin modes B^ K*^p'^ 
and (j)K* is of interest since the ~ 0.5 transverse polarization fraction observed in (pK* decays 
could also be due to contributions from new physics in the penguin loop. 

Table H summarizes the theoretical predictions from the naive-factorization model compared 
to the existing measurements of branching fractions and polarizations for the K*p and pp 
decay modes. The theoretical predictions for the branching fractions are taken from 6j. 
They have been updated in ^3] . The ranges given are obtained by varying the form factors 
and other parameters entering the calculation. The polarizations have also been predicted. 

Table 1: Comparison of the predictions from the naive-factorization model to the measure- 
ments for the branching fractions and polarizations in the K* p and pp decay modes. The 
branching fractions are in units of 10~^. The results for B^ —>■ K*^p'^ are those reported 
here. 



Mode B (10^^) B (10"^) B (10^^) Polarization Polarization 

prediction prediction [H| Measurement prediction j21 measurement 





6.6 


6 - 


10 


10.6;^« 


0.90 


0.96-:- 


K*+p- 


7.0 


6- 


10 


< 17.2 @ 90% CL 


0.90 




K*'^p+ 


9.0 


8 - 


12 


i7.o;|^ 


0.90 


0.79 ±0.09 




6.1 


7- 


12 


26.4111 


0.92 


0.89 ±0.07 


p+p 


24.0 


20 - 


-25 


27 ±9 


0.92 


98+° °2 

'-'•^"-0.09 



1.2 Angular analysis 

The analysis is done in the helicity frame (Fig. ^ as a function of the "helicity angles" 9k*o 
and 6p+. The angle 6k*o {Op+) is defined as the angle between the direction of the K*^ (p"*") 
and the direction of the tt^ (tt'^) coming from its decay in the vector meson's rest frame. An 
integration is performed over the angle between the vector-meson decay planes to simplify 
our analysis; this step is straightforward because the detector acceptance is independent of 
this angle. The longitudinal-polarization fraction can be extracted from the differential 
decay rate, parametrized as a function of 6k*o and 6p+ |15j : 

1 d^r 1 

^1 a J ^ " f^) sin^ Ok-o sin^ 6^^+ + /l cos^ Ok-o cos^ 6^^+. (1) 

1 a cos a K*od cos Up+ 4 



Experimentally, it is important to measure the branching fraction and the polarization simul- 
taneously because of their large correlation. The two decay products of the vector mesons 
have in general comparable momenta when transversely polarized, and asymmetric momenta. 
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Figure 1: Helicity frames for the vector- vector —>■ K decay. 



with one high- and one low- momentum decay product, when longitudinally polarized. As soft 
particles have lower reconstruction efficiency, the efficiency for reconstructing longitudinally 
polarized decays is about half that for transversely polarized decays. 

2 THE BABAR DETECTOR AND DATASET 

The results presented in this paper are based on data collected in 1999-2002 with the BABAR 
detector at the PEP-II asymmetric e~^e~ collider of the Stanford Linear Accelerator 
Center. An integrated luminosity of 81.85 fb~^ was recorded at the T(AS) resonance, cor- 
responding to 88.84 ± 0.98 million BB pairs. An additional 9.58 fb~^ data sample taken 
40 MeV below the T(4S') resonance is used in order to study the continuum background 
e~^e^ ^ qq {q = u, d, s, c). 

Charged particles are detected and their momenta measured with the combination of a silicon 
vertex tracker with five layers of double-sided detectors and a 40-layer central drift chamber, 
both operating in the 1.5-T magnetic field of a solenoid. Charged-particle identification 
(PID) is provided by measurement of the average energy loss (dE/dx) in the tracking devices 
and by an internally refiecting ring-imaging Cherenkov detector covering the central region. 
A K/ir separation of better than four standard deviations (a) is achieved for momenta 
below 3 GeV/c, decreasing to 2.5 a at the highest momenta found in B decay. Photons and 
electrons are detected by a CsI(Tl) electromagnetic calorimeter (EMC). 

3 ANALYSIS METHOD 

We reconstruct B^ K*^p^ candidates through the decays K*^ —>■ K^tt^ and p+ — > tt+tt'^, 
-jjO Inclusion of the charge conjugate processes is implied in this paper. 

The main backgrounds are from continuum events and BB events {B background). Event 
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selection reduces these backgrounds but as a maximum likelihood fit will be performed 
the selection is kept loose. Some non-resonant charmless B decays into four bodies are a 
particular problem as they have similar final states to our signal and their branching fractions 
are poorly known. 

3.1 Selection and discriminating variables 

Monte Carlo (MC) simulations ^\ of the signal and the B background as well as the off- 
resonance beam data were used to determine the selection criteria before the examination 
of the on-resonance beam data. 

The final state of the signal decay is K+tt+tt^tt^. We first select charged kaons and pions 
from charged tracks. Charged tracks candidates are required to originate from the interaction 
point: distance of closest approach to the interaction point less than 10 cm along the beam 
directions, and less than 1.5 cm in the plane transverse to the beam directions. We require 
that the charged pion candidates from the K*^ and decays not be identified as electrons, 
kaons, or protons, and that the charged kaon candidate from the K*^ decay agrees with a 
kaon hypothesis and be inconsistent with the electron and proton hypotheses. 

Next, we reconstruct tt^ candidates from photon pairs, where each photon has an energy 
larger than 50 MeV and exhibits a lateral profile of energy deposition in the EMC consistent 
with an electromagnetic shower ^HI- The n'^ candidate mass must satisfy 0.11 < m^^ < 0.16 



We then reconstruct K*^ and p"*" candidates. The mass of the and the p"*" candidates 
must satisfy \mK+n- - 0.896| < 0.125 GeV/c^ and |m^+^o - 0.769| < 0.375 GeV/c^. These 
mass windows correspond to about 2.5 times the resonance full width (Fig. IH}. Combi- 
natorial backgrounds dominate near the helicity angle regions |cos^i^*O p+| = 1, where the 
vector mesons tend to decay in a way that produces a soft particle. The effect is more 
important when this soft particle is a vr*^. If one assumes that the longitudinal polarization 
is large, as suggested by recent measurements of charmless vector-vector modes and theo- 
retical predictions, it is important to maintain the largest possible coverage in cos9k*o,p+- 
The requirements —0.95 < cos6'x*o < 1-0 (—0.8 < cos6'p+ < 0.95) are applied to reject 
candidates with a soft tt^ (tt^ and 7r+) coming from the K*^ (p+) decay; they allow the 
suppression of most of the combinatorial backgrounds while maintaining adequate efficiency 
for longitudinally-polarized decays. 

We select B candidates from the K*^p^ combinations using two nearly independent kinematic 



observables the beam energy-substituted B mass niEs = y (f + ViVbY / E"^ ~ Pb 
the energy difference AE = {EiEB — Pi-Ps — ^)/^/s, where ^/s is the beam energy in 
the T(4S') CM, and {Eb,Pb) and {Ei,pi) are the four-momenta of the B candidate and 
the e~^e~ initial state, both defined in the laboratory frame. For signal events, the rriEs 
distribution peaks at the S-meson mass and the AE distribution peaks near zero. We 
require B candidates to satisfy 5.21 < rriEs < 5.29 GeV/c^ and \AE\ < 0.15 GeV. When 
multiple B candidates exist in the same event, we select the one whose reconstructed vr*^ 



GeV/c2. 
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mass is nearest to the known vr*^ mass; we choose the candidate randomly from those that 
share the same 7r°. 

To discriminate signal from continuum background, we also use a neural network (NN) 
combining six variables: a Fisher discriminant made from two event-shape variables (see 
[T^): the cosine of the angle between the direction of the B and the collision axis (z) in the 
CM frame; the cosine of the angle between the 5-thrust axis and the z axis; the cosine of 
the angle between the 5-thrust axis and the thrust of the particles of the rest of the event; 
the angle between the direction of the 7r° and that of one of its daughter photons in the vr" 
rest frame (vr^ decay angle); and the sum of transverse momenta relative to the z-axis of the 
particles in the rest of the event. 

3.2 Sample composition 

The expected numbers of signal and background events in the data sample are summarized 
in Table 13 The B background is divided into categories, described below, that are modeled 
separately in the maximum likelihood fit (Sec. 13. 3p . 

Table 2: Number of events expected in the K*^p^ analysis for an 81.85 fh^^ data sample. 
A branching fraction of 20 x 10^^ and a longitudinal polarization of 75% are assumed for 
the signal. 



Category 


Expected number 




of events 


Signal 


167 


Continuum 


12400 


h —>■ charm backgrounds 


B+ ^ Z}07r+ 


409 


B+ ^ 5 V 


696 


B^ —>■ c, excluding D^n^ and D^p^ 


915 




904 


b —>■ charmless backgrounds 


p+p" (100% longitudinal polarization) 


3.2 




4.3 


"Five bodies" 


68 


"Other Charmless" 


291 



3.2.1 Signal 

The signal consists of correctly reconstructed "true signal" events as well as badly recon- 
structed "Self-Cross-Feed" (SxF) events. SxF events contain useful information as one of 
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the two vector mesons is usually correctly reconstructed. They are treated separately in the 
fit. 

The SxF events are due mostly to the misassignment of a charged (neutral) pion in about 60% 
(40%) of the cases. The SxF events are decomposed into three categories: events in which 
the K*^ is reconstructed correctly but not the p"*", events in which the p"*" is reconstructed 
correctly but not the K*^, and events in which neither vector meson is reconstructed correctly. 

The total signal selection efficiency, including both "true signal" and SxF events, is 12.6 ± 
0.1% (20.6 ± 0.2%) for longitudinally (transversely) polarized events. The fraction of SxF 
events is 25.1 ± 0.2% (10.9 ± 0.3%) for longitudinally (transversely) polarized events. 

3.2.2 B background 

Background from B decays can be split into that from b ^ c transitions and that from b 
charmless transitions. We further split the b ^ c background into four subcategories, each 
with its own term in the likelihood fit (Sec. 13. 3|) : B~^ D^n'^, B^ D^P^ ■, charm 
other than the first two, and B^ —>■ charm. The first two peak in mEs while the last two do 
not. In particular, B+ D°7r+, D° K+tx-'k^ and B+ D^p+, K+n- share the 

same final state K^n^-ri^ii^ with signal. 

The other charmless B backgrounds are modeled in four categories in the likelihood fit. The 
two B^ P^P^ and B^ —>■ K*^7i^ specific modes are separated due to their similarity with 
the signal. The rest of the charmless are divided into two additional categories. The "five 
bodies" category includes decay modes that have five particles in the final state and involve 
an intermediate resonance [a^,af ,aQ,aQ ,uj, fo, p'^ , p^, K*~^ , K*^] or one of the particles K~^, 
vr"'" or 77°. The "Other Charmless" category consists of all other charmless modes after all 
the above modes have been excluded. The modes in the "five bodies" category resemble 
the signal more than the ones in the "Other Charmless" category due to the resonances. 
Therefore their mEs distribution peaks slightly more at the 5-meson mass. However the 
yield of the "Other Charmless" category is four times larger (Table Ej). Since it is the most 
poorly known, it will be floated in the likelihood fit. 

None of the non-resonant charmless modes K*^tt^7t^, K^p^n^, and K^tt^tt^'k^ , with the 
same final state as the signal, have been studied experimentally. No theoretical model is 
currently available to predict their branching ratios and their decay kinematics. Therefore 
they are not modeled in the likelihood fit in this preliminary analysis. Yet they could produce 
substantial false signals since their final state is the same as for signal events. A systematic 
error will be assigned for this fact (Sec. 14. 4j) . 

To summarize, we have a total of 13 categories of events: the "true signal", three SxF 
categories, the continuum, four categories of charm i?-background, and four categories of 
charmless 5-background. In the fit, we float the sum of the yields in the "true signal" 
and the SxF categories, as well as the yields in the continuum and the "Other Charmless" 
categories. The yields in all the other 5-background categories are fixed. 



15 



3.3 Likelihood fit 



We perform an unbinned, extended maximum likelihood fit on the selected data sample 
to extract the signal yield, the longitudinal polarization f^, and the direct-CP- violating 
asymmetry Acp- Seven observables are used: tties, ^E, the helicity angles cos(6'i^.o) and 
cos(6'p+), the reconstructed masses of the two vectors mK*o and mp+, and the neural network 
output. 

All significant correlations are taken into account in the likelihood function described below 
or are covered by systematic errors (Sec. 14.2.11) . 

The extended likelihood function is then given by: 

A. / \ 



1=1 



^sigif^Sig' Psig,i + fsxFkPsxFk,i) + ^contPcont,i + 

y 5xFfc=l,3 B backg. J 



(2) 



where Ugig and ricont are the numbers of signal (including SxF) and continuum events that 
are floated in the fit. The numbers of events Uj in the i?-background category j are all fixed 
to their Monte-Carlo expectations, except in the case of the "Other Charmless" category for 
which the yield is floated. N'{N) are the expected (observed) total numbers of events in the 
data sample and i is the event index. Z^™^ and fsxFi^z are the fractions of "true signal" and 
of the three different SxF categories normalized to the total number of signal events nsig. 
These fractions are different in the longitudinal and transverse components and are taken 
from the simulation. 

The normalized P Probability Density Function (PDF) in each category is the product of 
the normalized PDFs of each observable, except for the continuum for which a joint PDF of 
the correlated helicities and vector-meson masses is used. 

PDF Models 

The irtESi Ai?, NN, cos^^x*", cos6p+, m^^.o, and mp+ PDFs of the "true signal" and the B 
background are taken from the simulation, but the means and widths of the signal Gaussian 
PDFs for tues and AE are corrected to account for the differences between data and Monte- 
Carlo observed in a D^tt^ (with D'^ — > K^tt^tt^) control samples. 

For the signal including the three SxF categories, the distributions of the vector masses 
and helicities are modeled using for each vector meson the distributions of equation (^Q) 
modified by a function Oi^-.o p+(^x*o,p+) proportional to the probability for a signal event to 
be reconstructed, correctly or not, and multiplied by the probability Pk*o p+i^K*o,p+) (1 — 
Pk*o p+{Ok*o^p+)) to reconstruct correctly (or not) the K*'^ or the p"*"; both functions depend 

on the vector-meson helicity angle. Different mass distributions -Pwi^1o^p+"^*'^(mx*o_p+) of the 
vector- meson mass mx*o^p+ are used whether the reconstruction of the meson is correct or 
not. For example, for the "true signal" longitudinal component we have: aK*o{9K*o)ap+{9p+) 
X Ppo{0K*o)Pl7^{ep+) X cos'^{9K*o)cos'^{ep+) X Pm]^,T{rnK*^)Pn{^^l''\mp+) . Another ex- 
ample is the transverse component of the SxF category with a correctly reconstructed 
K*^ described by: ai^.o(ei^.o)ap+(^p+) x P^™o'(^i^.o)(l - P*™'^(^p+)) x sin^ {6 k*o) srn^ {6 p+) x 
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Table 3: Summary of the multiplicative systematic uncertainties AB on the branching frac- 
tion B associated with the signal reconstruction efficiency and number of BB. 



Source AB 



TT^reconstruction 


8.6 


% 


Track reconstruction 


3.9 


% 


PID 


1.1 


% 


Number of BB pairs 


1.1 


% 


Total 


9.6 


% 



Pm^^o'^{mK*o)Pm^+''^^{mp+) . The model supposes that the reconstruction of the two vector 
mesons is independent. A systematic uncertainty for this assumption is included in the total 
uncertainty (Sec. I4.2.1|) . 

In the continuum, the mEs distribution is parametrized by an ARGUS function ^H], the /S.E 
distribution by a second order polynomial, and the neural network output (NN) distribution 
by the three-parameter function (1 — NNY""^^ a2(i-NN)+ a^ii-NN) )_ rjj^g shape parameters 
entering the PDFs of the ruEs, AE, and NN observables are also floated in the flt: the 
statistics of the "on-resonance" data is much larger than the "off-resonance" data, which 
is used for cross-checks. To model the correlation between the mass and the helicity of 
each vector meson, non-parametric PDFs are made from two 2-dimensional mass-helicity 
histograms extracted from the rriEs side-band {niEs < 5.25 GeV/c^ and A^A^ < 0.4), after 
subtraction of the remaining B background (less than 10%). 

4 SYSTEMATIC UNCERTAINTIES 

The systematic uncertainties are summarized in Tables El and EJ Table El displays the un- 
certainties on the efficiency of signal reconstruction and on the total number of BB pairs in 
the data set, each which contributes to a multiplicative uncertainty on the final branching 
fraction. The quadratic sum of the uncertainties in Table El is reported in Table 01 as "Signal 
reconstruction efficiency, number of BB pairs." It is one of the largest contributions to the 
overall systematic uncertainty on the branching fraction measurement. 

4.1 Fractions of "True signal" and SxF categories, and efficiencies 

Table m also gives the systematic uncertainties on the polarization fi, due to the uncertainty 
on the ratio of the selection efficiencies for the longitudinal and transverse polarization 
components (used to compute the effective polarization of the selected signal events). 

Other systematics on the signal yield and the polarization are due to the uncertainties on 
the relative fractions Z^™^ and fsxFi,2,3 of the "pure signal" and of the three SxF categories. 
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Table 4: Summary of the systematic uncertainties on the signal yield, on the branching frac- 
tion, and on the polarization. The asymmetric uncertainty from the non-resonant charmless 
backgrounds is presented separately from the other systematics. 



Source 


Ns B 


h 


Signal reconstruction efficiency, number of BB pairs 
Ratio of efficiencies for long. /trans, polarizations 
"True signal" and Self-cross-feed fractions 
PDF shapes 

NN shape in off-resonance data 
No. B backgrounds 


9.6% 

±0.0 ±0.0% 
±1.7 ±1.2% 
±6.4 ±4.5% 
±4.9 ±3.5% 
±1.8 ±1.3% 


±0.006 
±0.001 
±0.019 
±0.027 
±0.010 


Total 


±8.4 ±11.3% 


±0.035 


Non-resonant charmless backgrounds 
{K^n^Ti^Ti^ final state) 


-15.7 -11.1% 


±0.020 



which are taken from the simulation. Due to correlations between the K* and p detection 
and reconstruction, the fraction of events in one category of signal or SxF is slightly different 
from the product of the probabilities Pj^^o p+{dK*o^ p+) (Sec. 13. 3|) to reconstruct correctly or 
not the vector meson, averaged over the helicity angles. The difference between the results 
obtained with this set of signal and SxF fractions and the set of values taken directly from 
the simulation is taken as the systematic uncertainty. 

4.2 Uncertainties on PDFs shapes 

The systematic labeled "PDF shapes" in Table|3]is due to the uncertainty of the PDF shapes 
for the signal, continuum, and B backgrounds. It includes the error on the correction of the 
mean and width of the "true signal" and AE gaussian distributions for the differences 
observed between the data and the simulation (Sec. 13. 3|) . It includes also systematics 
associated with corrections of biases from the PDF models for the signal, SxF, and continuum 
(Sec. I4.2.1|) . and uncertainties due to the limited statistics in determining PDF shapes (Sec. 

4.2.1 Systematics from the PDF model itself 

We first checked that the cut cos^p+ > —0.8 is sufficiently tight not to leave low momentum 
7r° backgrounds that are not properly simulated: the results are stable when the analysis is 
repeated after a tighter selection on the helicity angle (cos6'p+ > —0.5). 

The model for the "true signal" and the SxF is tested using fully simulated signal, divided 
into sub-samples embedded in toy simulations of the other backgrounds. Each sample has 
the same number of signal events and the same polarization as is fitted in data. A small 
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bias of 4.1 ± 1.4 events is observed on the signal yield as well as a (—1.2 ± 0.5)% bias on the 
polarization. 

In the continuum, the masses and helicity distributions are modeled by using 2-dimensional 
mass-helicity histograms extracted from the mEs side-band (Sec. IH.Hj) . A slight bias of 
—2.2 ± 0.2 events on the signal yield and of (—0.9 ± 0.1)% on the polarization is associated 
with this model. This is estimated using toy Monte-Carlo samples, generated with the PDF 
used for the data fit. For each toy experiment new PDFs of the masses and helicities of 
the vector mesons are constructed following the same procedure as in the real data sample. 
The results obtained when using the new PDFs are slightly biased compared to the results 
obtained when using the PDFs the toy sample was generated from. 

The two biases from the signal and continuum models are corrected for in the final result, 
and conservative systematic errors, each equal to half of the correction, are assigned. 

4.2.2 Systematics from the limited statistics to determine the PDFs 

The systematic uncertainty associated to the shape of the parametrized PDFs is determined 
by varying the parameters within their statistical error obtained from a fit on the full sim- 
ulation. The non-parametric PDFs are varied by generating toy samples using the original 
PDFs used in the real data fit. For each toy sample, new PDFs are made out of the generated 
distributions of the observables. The systematic error is the dispersion of the differences in 
the fit results between the new and the original PDFs. 

4.2.3 Neural network shape parameters in the continuum 

For the continuum, the itles, ^E, and NN PDFs are parametrized. The parameters are 
determined from real data (Sec. 13. 3|) . Their values agree with the values fitted on the 
off-resonance data for the mEs and AE distributions, but not for the NN. A systematic 
uncertainty is assigned, labeled "NN shape in off-resonance data" in Table |3 It is equal 
to the difference on the results when using the two sets of parameters. Though it is an 
uncertainty on a PDF shape, it is displayed separately from the "PDF shapes" category as 
it is relatively large. 

4.3 Number of events in the 5-background categories 

Other systematic uncertainties, labeled "No. B backgrounds" in Table EJ come from the 
uncertainty on the numbers of events in each 5-background category that are fixed in the 
fit. The branching ratio of the —* p~^p^ decay is varied within the errors of the combined 
measurements of [3' and 0: 3(3^^ p^)= (26.4^54) x 10~^. The unknown branching 

fraction of the i^*°7r° background mode is varied by 100%, assuming a central value of 
7.5 X 10~^. This value is half of the value of the measured branching fraction for the K*^tt^ 
decay mode pO (f^i 15 x 10~^), based on the isospin rule and the assumption of penguin 
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dominance. The number of events is varied by 20% in all the b ^ c categories and by 50% 
in the "Five bodies" category of charmless- 1? background. Note that there is no systematic 
error assigned for the "Other Charmless" category of events as this yield is floated in the fit 
to data. 

4.4 Contribution from non-resonant charmless backgrounds de- 
caying into A'+TT+TT TT^ 

The non-resonant charmless 5-decay modes K*^tt^tt^, K^p'^tt^, and K^tc^tt^tt^ have the 
same final state as the signal and are not modeled in the fit. The associated systematic 
error is estimated by the difference in the data fit result when the yields of these background 
modes are floated or fixed to zero: PDFs of these modes can be constructed from their full 
simulation, based on a simplified phase space model, but the available statistics is limited 
to a few hundreds of events. We obtain an asymmetric error of —15.7 events (—11.1%) 
on the signal yield, which is systematically overestimated when these background modes 
are not modeled in the fit. We get an additional symmetric systematic error of 2.0% on 
the polarization. This main systematic error is a preliminary estimation and is presented 
separately from the other systematics, with the label "non-resonant" . 

Note that this study seems to show that most of these backgrounds are negligible, probably 
because their decay products populate different regions of the phase space. The only channel 
which has overlapping phase space is {K~^'7i~)s-waveP~^- A more detailed study of these 
backgrounds will be necessary. 

5 PHYSICS RESULTS 

We obtain from the fit 147.3l22'3(stat) signal events and a longitudinal polarization frac- 
tion fi = 0.77 ± 0.08(stat). We then correct for the two small biases described in Sec- 
tion l4.2.r| from the models of the signal and self-cross-feed, and of the continuum. We obtain 
141.0;!;22!3(sfat) ± 15.9(syst)^°5°7(non-resonant) signal events and Jl = 0.79 ± 0.08(stat) ± 
0.04(syst) ± 0.02(non-resonant). The polarization observed is consistent with both the po- 
larization found in the other pure penguin modes such as K*(f) {fi ~ 0.5) and with purely 
longitudinal polarization (/^ ~ 1). From the number of signal events, the fraction of longi- 
tudinal polarization, the selection efficiencies determined for the transverse and longitudinal 
polarization components, and the branching fractions B{K*^ — > K^tt^), i3(7r° — > 77), we 
compute the branching fraction: 

K*^p+) = [17.0 ± 2.9(stat) ± 2.0(syst)l°:° (non-resonant)] x 10"^- 

The impact of the uncertainties on B{K*'^ K^tc^) and B{tt^ — > 77) is negligible compared 
to the other systematic errors. The statistical error on the branching fraction results from 
the statistical errors on the signal yield and the polarization, taking into account their 
correlation (Fig. |2I). The systematic uncertainty on the branching fraction results from the 



20 




1 



5 10 15 20 25 30 35 
Branching Fraction (x 10 ) 



Figure 2: Countours at the 1, 2, 3, 4, and 5 cr levels showing the correlated statistical 
uncertainty on the branching fraction and on the polarization. 



propagation of the systematics uncertainties on the signal yield and on the polarization. The 
systematic error on the signal yield (resp. polarization) of 15.9 events (resp. 4%), results in 
an systematic error of 1.9 x 10~^ (resp. 0.4 x 10~^) on the branching fraction. Finally, the 
"non-resonant" systematics on the signal yield is propagated to the branching fraction. 

This result is consistent with the isospin rule predicting that B{B'^ —>■ K*^p^) ^ 2x B{B^ — > 
K*~^p^), assuming the dominance of gluonic penguin diagrams in the K*^p^ decay 

and the measured B^ — > K*^p'^ branching fraction of [10. 61^35] x 10~^ This result is 
about 2 standard deviations away from the predictions from naive factorization models jHlE] 
(Tabled). 



It is consistent with zero, as expected for this pure penguin mode. The systematic uncertain- 
ties on Acp arises from a limit on the possible size of charge- dependent tracking efficiency 
and particle identification biases. The uncertainty from the tracking is the linear sum of the 
0.3% uncertainty for each charged track. The uncertainty from the particle identification is 
estimated at 4%. 

The distributions of rriEs, ^E, and the vector-meson helicities and reconstructed masses are 
shown in FigureElwith the signal enriched by selecting events with large signal to background 
likelihood ratios on the discriminating variables not shown in each plot. 



We also extract the direct-CP-violating charge asymmetry: 



Acp{B 



K*^p+) = [-14 ± 17(stat) ± 4(syst)]%. 
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6 SUMMARY 



We measure the branching fraction and the fraction of longitudinal component for the decay 
— * K*^p^ using a maximum likelihood technique. We use a data set corresponding to a 
total integrated luminosity of 81.85 fb~^ taken on the T^iS) peak. A signal is observed for the 
first time with a significance of greater than 5a. Prom a fitted signal yield of 141.0l22!3(stat)± 
8.4(syst)_j^5 ^(non-resonant) events we obtain the branching fraction: 

B{B+ K*^p+) = [17.0 ± 2.9(stat) ± 2.0(syst)+$:^(non-resonant)] x 10"^ 

the longitudinal polarization fraction fi' 

Jl = 0.79 ± 0.08(stat) ± 0.04(syst) ± 0.02(non-resonant), 

and the direct- CP- violating asymmetry: 

Acp{B^ K*^P'^) = [-14 ± 17(stat) ± 4(syst)]%. 
These results are preliminary. 
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Figure 3: Distributions of rriEs, Ai?, K*^ and reconstructed masses and helicities in the 
B'^ K*^p'^ decay, enhanced in signal component by selecting events with large signal to 
background likelihood ratios on the discriminating variables not shown in each plot. In the 
top four plots, the black dots are the data, the blue plain line is the fitted distribution for the 
full data sample, the black dashed line is the fitted distribution for all backgrounds, and the 
dotted red line is the distribution of only the continuum events. The two plots in the bottom 
show the distributions of the helicity angles. The data (black dots), after subtraction of the 
backgrounds, are compared to the signal component of the fitted PDF (plain black curve). 
In dashed lines are the results expected for different fractions of longitudinal component: 
100% (blue), 50% (red), 0% (cyan). The pure transverse case is excluded by the data. 



23 



References 



[1] A. Buras, R. Fleischer, S. Recksiegel, and F. Schwab. Anatomy of prominent B and K 
decays and signatures of CP-violating new physics in the electroweak penguin sector. 
(2004). hep-ph/0402112. 

[2] M. Neubert and J.L. Rosner. Determination of the weak phase 7 from rate measurements 
in ttK^tttt decays. Phys. Rev. Lett, 81:5076-5079, (1998). hep-ph/9809311. 

[3] B. Aubert et al. [BABAR Collaboration]. Rates, polarizations, and asymmetries in 
charmless vector-vector B meson decays. Phys. Rev. Lett., 91:171802, (2003). hep- 
ex/0307026. 

[4] J. Zhang et al. [BELLE Collaboration]. Observation of B^ P^P^- Phys. Rev. Lett., 
91:221801, (2003). hep-ex/0306007. 

[5] B. Aubert et al. [BABAR Collaboration]. Search for the decay B^ K*+ p' . (2004). 
BABAR-CONF-04/041, SLAC-PUB-10636, hep-ex/0408035. 

[6] A. Ah, G. Kramer, and Cai-Dian Lu. Experimental tests of factorization in charmless 
non-leptonic two-body S decays. Phys. Rev., D58:094009, (1998). hep-ph/9804363. 

[7] R. Aleksan, P. F. Giraud, V. Morenas, O. Pene, and A. S. Safir. Testing QCD factor- 
ization and charming penguins in charmless B — > PV. Phys. Rev., D67:094019, (2003). 
hep-ph/0301165. 

[8] A. Kagan. Polarization va. B ^ VV decays. (2004). hep-ph/0405134, submitted to 
Phys. Lett. B. 

[9] B. Aubert et al. [BABAR Collaboration]. Observation of the decay B^ — > p~ and mea- 
surement of the branching fraction and polarization. Phys. Rev., D69:031102, (2004). 
hep-ex/0311017. 

[10] B. Aubert et al. [BABAR CoUaboration] . Measurement of the B^ (t)K*^ decay 
amplitude. (2004). BABAR-PUB-04/031. 

[11] A. Bosek et al. [BELLE Collaboration]. Phys. Rev. Lett, 91:201801, (2003). hep- 
ex/0307014. 

[12] Y. Grossman. Int J. Mod. Phys., A19:907, (2004). 

[13] P. Colangelo, F. De Fazio, and T. N. Pham. The riddle of polarization in S — > VV 
transitions. (2004). hep-ph/0406162. 

[14] G. Schott. Study of the rare decays B K* p/ pp and search for CP violation in these 
modes in the BABAR experiment. Ph.D. Thesis, Paris 6 University, July (2004). 

[15] G. Kramer and W.F. Palmer. Phys. Rev., D45:193, (1992). 

[16] B. Aubert et al. [BABAR Collaboration]. Nucl. Instrum. Methods, A479:l-116, (2002). 



24 



[17] S. Agostinelli et al. [GEANT4 Collaboration]. Nucl. lustrum. Methods, A506:250, 
(2003). 

[18] B. Aubert et al. [BABAR Collaboration]. Phys. Rev. Lett, 89:281802, (2002). 

[19] H. Albrecht et al. [ARGUS Collaboration]. Z. Phys., C48:543, (1990). 

[20] B. Aubert et al. [BABAR Collaboration]. Measurements of the branching fractions of 
charged B decays to K^n^n^ final states. (2003). hep-ex/0308065, Submitted to Phys. 
Rev. D. 



25 



